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MICROLENSING EVIDENCE THAT A TYPE 1 QUASAR IS VIEWED FACE ON 
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ABSTRACT 

Using a microlensing analysis of 11-years of OGLE V-band photometry of the four image gravi- 
tational lens Q2237+0305, we measure the inclination i of the accretion disk to be cosi > 0.66 at 
68% confidence. Very edge on (cos? < 0.39) solutions are ruled out at 95% confidence. We measure 
the V-band radius of the accretion disk, defined by the radius where the temperature matches the 
monitoring band photon emission, to be Ry = 5. x 10 15 cm assuming a simple thin disk model 
and including the uncertainties in its inclination. The projected radiating area of the disk remains 
too large to be consistent with the observed flux for a T oc i?~ 3 / 4 thin disk temperature profile. There 
is no strong correlation between the direction of motion (peculiar velocity) of the lens galaxy and the 
orientation of the disk. 

Subject headings: gravitational lensing — methods: numerical — quasars: general — quasars: indi- 
vidual (Q2237+0305) - 



1. INTRODUCTION 

In the AGN unificat ion model (e.g. lAntonuccl 119931 : 
lUrrv fc Padovanil [1995) , orientation differences among 
intrinsically similar objects are thought to account for 
many of the different observational properties of AGN. 
In particular, a dusty "torus" may frequently obscure the 
central engine from direct observation when viewed edge- 
on. The bright, Type 1 broad line quasars are thought 
to be viewed mostly face-on, so that the accretion disk 
is not obscured by material in the equatorial plane, and 
Type 2 narrow line quasars are viewed closer to edge-on 
through the obscuring material. Thus, if we could resolve 
the disk of a Type 1 quasar, we would expect it to be 
closer to face-on than edge-on. Unfortunately, familiar 
methods cannot resolve accretion disks, so we have only 
indirect measures of AGN disk ori entation. For exam- 
ple, t he projected axes of r adio jets dBlandford fc Konigll 
I1979D and ionization cones (|Elvisll2000l ) both support this 
picture. There are, however, no actual measurements of 
disk orientation. 

While quasar accretion disks are too small to be re- 
solved by direct imaging, gravitational microlensing pro- 
vides a natural telescope to study the structure of quasar 
accretion disks and the properties of cosmologically dis- 
tant lens galaxies w here we see multiple images of back- 
ground quasars fsee iWambsganssI 120061) . In addition to 
the mean potential of the lens galaxy, each image is also 
magnified by the microlensing effects of the nearby stars. 
Since the observer, the lens galaxy and its stars, and the 
quasar are all moving, microlensing is observed as un- 
corrected time variability in each of the quasar images. 
The amplitudes of these variations depend on the struc- 
ture of the accretion disk and the properties of the lens 
galaxy. 

Quasar microlensing is most sensitive to the projected 
area of the accretion disk relative to the source plane 
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where G is the gravitational constant, c is the speed 
of light, (M) is the mean stellar mass of the stars, 
-Dls>-C>OL, and Dos are the angular diameter distances 
between the lens-source, observer-lens, and observer- 
source respectively, and we have used the lens and source 
redshifts for Q223 7+0305 (zi = 0.0394, z s = 1.685, 
iHuchra et ail 119851 Q2237 hereafter). The smaller the 
accretion disk, the higher the variability amplitude from 
microlensing. In general, the emission profile is difficult 
to determine, as models having simil ar half light radii 
show similar microlensing variability (jMortonson et all 
2005). There has been little examination of other struc- 
tural parameters of disks, except for iCongdon et al.l 
(|2007|) . who demonstrated in simulations that the mi- 
crolensing signal is sensitive to the ellipticity and orien- 
tation of the accretion disk. If detectable in practice, 
measuring the apparent ellipticities of accretion disks 
provides an important test of AGN unification models 
and opens the possibility of examining the complex ef- 
fects of relati vity on the apparen t surface brightness of 
the disk fe.g. lHubenv et al.|[200ll) . 

Measurements of quasar disk sizes using microlensing 
are now c ommon. Recent effo rts have studied individual 
sizes fe.g. lMorgan et alj|2008l). the relationsh i ps between 
size and wavelength (lAnguita et all 120081: iBate et al.l 



2008: | Eigenbrod et al.l 



2008ri iPoindexter et al.l 120081: 
Floyd et all 120091: iMosauera et all 120091) . size and black 
hole mass (jMorgan et al.l 12010; , and the sizes of ther- 
mal and non-thermal emission regions (Poole y et al 



20071 iMorgan et al. 2008b; Chartas et all2009H Dai et a. 



d. 
d. 



20091) . All these studies used circular accretion disks 
and static magnification patterns that neglect the ran- 
dom motion of stars in the lens galaxy. 




Fig. 1. — Example of a trial source trajectory (dark line segments) superposed on instantaneous point-source magnification patterns for 
(M/Mq) = 0.3. Darker shades indicate higher magnification. An HST H-band image in the center labels the images and the corresponding 
magnification patterns. Each pattern is rotated to have the correct orientation relative to the lens. This particular LC2 trial has an effective 
lens- plane velocity of ~ 600 km s -1 Northeast. The "bow tie" (top) exhibits the two peaks in the major-axis position angle distribution 
together with the shaded region it approximately represents the 68% confidence region. The solid disk (right) has a radius of 10 17 cm. 



Determining the shape and orientation of a disk de- 
pends on the correlations between the anisotropic struc- 
ture of the disk and the anisotropic structure in the 
magnification patterns created by the sh ear (tidal grav- 
ity) l ocal to each image (see Figure Q] iCongdon et al.l 
2007). Existing microlensing studies cannot safely ex- 
plore these issues because they neglect the motions of 
the stars in the lens galaxy and use "static" magnifi- 
cation patterns. Since the stellar velocity dispersions 
of lens galaxies are comparable to the peculiar veloci- 
ties of galaxies, the patterns change on the same time 
scale as the source traverses them. Ignoring these stellar 
motions will overest imate the coherence of the magnifi- 
catio n patterns (see lWvithe et al.ll2000t iKochanek et al.l 
I2007J ) and likely render estimates of disk shapes unreal- 
istic. With a few exceptio ns that do no t focu s on disk 
structure (see Paper I and lWvithe et all ([1999D ). analy- 



ses of microlensing data have used static magnification 
patterns because of the computationa l challenges. In 
Paper I (Poind exter fc Koch anck 200jJ, we solved these 
computational problems and examined the peculiar ve- 
locity of the lens galaxy of Q2237 and the mean mass of 
its stars. In this paper, we measure the size, inclination, 
and position angle of the accretion disk in Q2237. In fJ2] 
we describe the data set used for our microlensing anal- 
ysis, our disk model, and outline our overall approach. 
Our results are presented in <J2 and a discussion follows 



2. DATA AND METHODS 

We analyze the nearly 11 years of Optical Gravitational 
Lensing Experiment ( OGLE) V-band pho tometric moni- 
toring data for Q2237 (jUdalski et al.ll2006[ ). To speed our 
analysis and as a cross check on the results, we divided 
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the data into two separate light curves. The first light 
curve (LCI hereafter) ranges from JD 2,450,663 to JD 
2,452,621 and consists of 100 epochs. The second light 
curve (LC2 hereafter) has 230 epochs from JD 2,452,763 
to JD 2,454,602. Each light curve covers just over 5 
years. We broaden the OGLE uncertainties by our esti- 
mate of the systematic uncertainties in the photometry 
of 0.02, 0.03, 0.04, and 0.05 m agnitudes in quadrature 
(jPoindexter k, Kochanekll2009fl . 

We analyze these lig ht curves using the Bayesian 
Monte Carlo method of iKochanekl (|2004ft . expanded to 
include motions of stars as detailed in Paper I. For each 
epoch of the light curve we generate a magnification pat- 
tern including the random motion of the stars. We used 
fixed mean masses of (M) = 0.01,0.03,0.1,0.3, 1,3, and 
10 M & , and a mass function of dN/dM oc l" 1,3 with a 
dynamic range M max /M min = 50 based on lGouldl §000). 
In Paper I we find that the best fit to the data is for 
(M) = 0.3A1©. The stars are assigned a random ve- 
locity in each coordinate from a Gaussian distribution 
with a = 170 km s" 1 based on the observed velocity 
dispe rsion (van de Ven , G. 2009, personal communica- 
tion, iTrott et al.l (|2008t )). The orbit of Earth (parallax 
effect) and the rotation of the lens galaxy are both in- 
cluded in the simulation. We convolve the patterns with 
the disk models described in $2j We draw the bulk ve- 
locities of the observer, lens galaxy, and source from 
a Gaussian of dispersion a = 1000 km s -1 on the lens 
plane in each coordinate. We later reweight the results 
to a more compact velocity prior based on the projec- 
tion of the CMB dipole velocity (jHinshaw et al.1 12009ft 
onto the lens plane, which is small (— 50, — 23) kms^ 1 , 
and the (ID) peculiar velocity dispersions of the lens 
and source of 327 km s -1 and 23 kms" 1 (estimated from 
iTinker. Wetzel. fe~Z ehavi 2009j), respectively. Because of 
the low lens redshift and the small projected dipole, the 
peculiar velocity of the lens is by far the most impor- 
tant factor (see Paper I). We then randomly draw light 
curves for each image and fit them to the data. Bayes 
theorem is used to combine the goodness of fit for the 
trials as measured by a \ 2 statistic, into probability dis- 
tributions for each variab le of interest. The se procedures 
are described in detail in IKochanekl (|2004l ) and Paper I. 

2.1. Accretion Disk Model 

We employ a generic thin disk model for which the 
sur face temperature scales as T cx i?~ 3 / 4 with radius 
R (|Shakura fc Sunvaevl 119731 ). The microlensing sig- 
nal is primarily sensitive to the half-light radius of the 
disk (Mortons on et al.ll2005D . which controls the effec- 
tive smoothing area of the disk, and, to date, studies 
have been unable to distinguish differing radial profiles 
(e.g.|Kochanck200l). A more realistic disk model would 
include a central hole whose size depends on the last 
stable orbit and general relativist ic effects modify th e 
underlying intensity profiles (e.g. iHubenv et al.l l200ll ). 
We have not presently pursued these features because 
they add many model parameters, and because it is un- 
clear whether the temperature profile of the thin disk 
model is correct. In particular, microlensing studies 
suggest the need for a flatter temperature profile (e.g. 
Pooley et aT1l2007t IMorgan et al.lfeOlOl; iPoindexter et al.l 
20081 Eigenbrod et al.l 12008b) in order to reconcile the 



microlensing sizes with the observed optical fluxes. Such 
changes in the temperature profile have also been sug- 
gested to explain the deviations of the observed spectral 
slope from the th in disk emission (see t he reviews by 
IKoratkar fc Blaesl (fl999ft and lBlaesl (|200l V 

Thus, we assume that the face-on (cosi = 1) surface 
brightness of the disk is 
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corresponds to the radius where the disk temperature 
equals the photon energy, kT — hc/X rcs t, Mbh is the 
black hole mass, M is the accretion rate, L/L-g is the 
luminosity relative to Eddington luminosity, and r/ = 
L/(Mc 2 ) is the radiative efficiency of the accretion disk. 
The half light radius of the disk is R 1/2 = 2.UR X . We 
now simply treat the disk as an infinitely thin disk viewed 
at inclination angles i, selected from a uniform distribu- 
tion in cosi = 0.2,0.4,0.6,0.8, and 1.0 (face-on). We 
include no relationship between the surface brightness 
and the viewing angle other than this simple projection 
effect in this first exploration of the problem. For full rel- 
ativistic disk models, there are many complexities such 
as Doppler shifts, red shifts, and bending of ray trajec- 
tories (e.g. I Agoll fl997t ) . but the bulk of the optical flux 
comes from relatively large radii, Ry/Rg ~ 30 for R g = 
GA/bh/c 2 = 2x 10 14 c m, where these effect s are less im- 
portant. Here we use [Morgan et all (|2010ft 's estimated 
black hole mass of 1.3 x lO 9 ]^ found by apply ing the 
virial relation of iVestergaard k, Peterson! (12006ft to the 
C IV line width measurement from lYee fc De Robertisl 
(1991). Since the orientation of the projected disk rela- 
tive to the magnification patterns also affects the results, 
we considered 18 major axis position angles from 0° to 
170° in steps of 10°, where the remaining angles are cov- 
ered by the reflection symmetry of the disk model. We 
parametrized the size of the disk by the projected area, 
7ri?yCOSz, as this is likely to have less correlation with 
the inclination angle because it keeps the projected area 
of the disk constant. 

We must also worry about whether all the observed 
emission arises directly from the accretion disk. The 
light curves can be contaminated by broad line emis- 
sion on much larger scales or some of the emission 
from the disk can be scattered on large r scales (see 
IDai et al.ll2009l : IMorgan et al.H2008bL l2T)10ft . At V- band, 
the contamination from broad line emission, mainly Fe 
pseudo-continuum emission and C III] A1909, contributes 
of order 20% of the fl ux in the spectral models of 
lEigenbrod et~aTl (|2008aft . To examine the effect of this 
dilution we ran models with 0%, 20%, and 40% contam- 
ination by emission on large, un-microlensed scales for 
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Projected Area [log(cm 2 )] 

Fig. 2. — The projected area 7r_RyCosi distribution of the ac- 
cretion disk. The blue dotted (red dashed) curve shows the re- 
sults from the analysis of only LCI (LC2). The black curve is 
the combined result from both LCI and LC2. The horizontal 
bars compare 68% confidence regions from earlier studies assuming 
(M) = 0.3Mq for estimat es which depended on the mean mass. 
The upper limit found by IWvithe et all COOObl) is at 99% con- 
fidence. The "flux" estimate is the area predicted based on the 
observed flux (Equation and is independent of inclination. The 
"theory" estimate is the area thin disk theory predicts assuming 
cos i = 1, based on the estimated black hole mass for Q2237 (Equa- 
tion 

the (M) = 0.3Mq case. Adding an unmicrolensed con- 
tamination fraction, /, in the V-band also decreases the 
flux size estimate by (1 — f) 1 ^ 2 unless the contamination 
is due to scattering o f the disk emission on large scales 
(see the discussion in M organ et al.|[2010f) . 

3. RESULTS 

We estimate the projected area of the accretion disk 
7ri?yCosi, the deprojected radius Ry, the disk inclina- 
tion i, and the major axis position angle. We also do 
a limited set of tests with different amounts of unmi- 
crolensed flux that may influence our analysis. We quote 
the results from the combined analysis of LCI and LC2, 
but also show the results from the independent analyses 
of LCI and LC2. Since the results are always mutually 
consistent, we only report quantitative results for the 
combined analysis. 

We find (Figure [5]) the projected V-band area of the 
quasar defined by 7ri?yCOS« to be 8.7 x 10 30 cm 2 < 
7ri? v cos i < 5.5 x 10 31 cm 2 (4.7 x 10 30 cm 2 < ttR\ cos i < 
9.4 x 10 31 cm 2 ) at 68% (95%) confidence with a median 
of 7T.RyCOsi = 2.2 x 10 31 cm 2 , where the size scale Ry 
is defined by Equation [3] After de-projecting the area 
based on each trial's inclination, we find that the V-band 
radius of the accretion disk is 3.5 x 10 15 cm < Ry < 
1.0 x 10 16 cm (2.6 x 10 15 cm < Ry < 1.4 x 10 16 cm) 
at 68% (95%) confidence with a median value of Ry = 
5.8 x 10 15 cm (Fi gure Et. This is consis tent with our 
earlier results in [Kochanck] ()2004j ) and iMorgan et al.l 
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Disk Scale Length [log(cm)] 

Fig. 3. — The disk scale length Ry after de-projection from the 
joint analysis of LCI and LC2. The blue dotted (red dashed) curve 
shows the results from the analysis of only LCI (LC2). The earlier 
results are shown assuming cosi = 1 for these face-on models. A 
scale on the right shows how they would shift to larger sizes if 
cos i < 1 . The flux size has the same inclination dependence as the 
earlier results, while the theory size is independent of inclination. 
Our test with 20% (40%) of light being contamination emitted on 
much larger scales resulted in a 26% (55%) smaller radius. 



(2010), of Ry = 3.7± 2 ;? x 10 15 / y/cosicm. and Ry = 

3.2l 3 g x 10 15 /Vcosicm using the same method without 
dynamic patterns, a smaller velo city prior, and shorte r 
light curv e s. Ot h er analyses by IWvithe et al.l (l20 00b) , 
lYoneharal (I2001T). IVakulik et all (I2007D . lAnguita et al.1 
(20081). and lEigenbrod et al.l (pOOSbTT have found gener- 
ally consistent results of Ry < 2.5 x 10 15 cm (at 99% 
confidence), Ry < 1.8 x 10 17 cm, Ry - 2 x 10 15 cm, 
Ry = 8.0±tl x 10 15 y/M/0.3M Q cm, and Ry = 3.3±|^ x 
10 15 cm (see Figure [3]), but using less data and with far 
stronger systematic assumptions. These estimates are all 
for cosi = 1 and should scale as 1/V cos i. If we compare 
projected areas (Figure [2]), then the comparisons are (to 
1st order) independent of the inclination angle, while if 
we compare disk scale lengths, there will be an inclina- 
tion angle dependence. While generally consistent with 
our results, our calculations use more data and are con- 
siderably more realistic, making it difficult to evaluate 
differences, especially the uncertainties. 

We can compare our measurement to those predicted 
by the observed fluxes or thin disk theory. If we as- 
sume only t he T oc R~ 3 / 4 temperatur e profile of a thin 
disk model (jShakura fc Sunvaevlll973T ) , where we can ig- 
nore the inner disk edge for these wavelengths, then the 
observed flux constrains the disk size by matching the 
integrated flux from a disk with the emission profile of 
Equation [3] to the observed flux. We estimate that the 
magnification corrected I-band flux is / = 18.03 ± 0.44 
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mag, and this corresponds to a disk radius of 
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where -Dos Ah is the angular diameter distance to 
the quasar relative to the Hubble radius, rn (see 
iMorgan et all 12010ft . This gives ttR\ cos 
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Assuming our best fit inclination, cos i 
_ x 10 15 cm. 

lAgol et ahl ()2009f ) estimate that the bolometric lumi 
nosity of the quasar is LAgol = 4 x 
corresponds to a Eddington factor of 



10 46 ergs s 1 , which 
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that is typical of luminous quasars ([Kolfmeier et al.l 
2006|). If we use this to replace the L/Le factor in Eqn. 
[3l the size estimate becomes 
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which now depends relatively weakly on the black hole 
mass. The uncertainties in these estimates are logarith- 
mic, corresponding to 0.4 and 0.3 dex respectively for 
Eqns. [3] and [5] if the masses, Eddington factors, lumi- 
nosities and efficiencies are viewed as being uncertain by 
a factor of 3. We find that the disk is large compared to 
the estimate based on the observed flux, although this 
is modestly reduced by the inclination, and small com- 
pared to a thin disk radiating close to Eddington with 
?/ = 10% efficiency. This is a discrep ancy common to all 
microlensing estima tes at present (see lPoolev et alJfeOOTt 
IMorgan et aU[20lol ). 

Where these size estimates are incremental improve- 
ments over earlier results from using additional data and 
an improved physical model, our results for the inclina- 
tion and orientation of the disk are entirely new. The 
preferred inclination is cosi > 0.66 at 68% confidence 
where cos % = 1 is face-on (Figure [4j . Such a relatively 
face-on inclination is consistent with the expectations of 
the AGN u nifica tion model. We can compare this to 
lAgol et ail (|2009ft 's model for the mid-infrared spectral 
energ y distribution of Q2237 u sing the du s t toru s mod- 
els of lFritz et all (|200l . While lAgol eTail ([2009ft do not 
trust the quantitative results, their preferred viewing an- 
gles (i < 19°, FigureHJ) are consistent with our inferences 
for the inclination angle of the disk. 

Figure [5] shows our estimate of the position angle (PA) 
of the major axis of the disk, and Figure Q] shows the 
68% confidence range of this PA for comparison to the 
anisotropies in the magnification patterns. We find con- 
sistent estimates of this position angle from both LCI 
and LC2 (Figure [5]). The preferred orientation is for 
the major axis to be roughly parallel to the shear in 
image B, but perpen dicular to it in images C and D. 
iCongdon et al.l ()2007ft found that microlensing variabil- 
ity is enhanced when the major axis is aligned with the 
shear, so our estimated alignment helps to explain the 
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Fig. 4. — Cosine of the disk inclination (cos i = 1.0 is face-on). 
The blue dotted (red dashed) curve shows the results from the 
analysis of only LC I (LC2). The p robability at cos i = was 
denned to be 0. The lAgol et al.l (120091 ') estimate is based on models 
of the mid-IR SED and, according to the authors, should not be 
interpreted quantitatively. 




50 100 150 

Disk Position Angle [degrees] 

Fig. 5. — Major axis position angle (North through East) of the 
accretion disk. The blue dotted (red dashed) curve shows the re- 
sults from the analysis of only LCI (LC2). The shear position 
angles for each image are indicated by the label arrows. The pre- 
ferred orientation of the major axis is to be parallel to the shear of 
images A and B but perpendicular to C and D. The face-on solu- 
tions are not included in these distributions because they provide 
no information on the position angle. 
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Fig. 6. — Major axis position angle of the disk relative to the 
direction of motion (N through E PAs). Here we exclude the face- 
on trials since they add no useful information to this distribution. 

higher variability of images A and B, and the lower vari- 
ability of D, but not the variability of image C. For this 
lens there are no other observations to which we can com- 
pare the orientation. However, if we had observations 
of the quasar host galaxy, we could compare the PA of 
the galaxy to that of the disk as a check on the rela- 
tive orientations of the angular momentum vec tors of the 
accret i on di sk and the galaxy. For example, lYoo et aT] 
( 20051, f2006h used images of lensed quasar hosts to con- 
strain the host axis ratio and position angle in 4 lensed 
quasars (to roughly < 20°). While both the peculiar 
velocity direction of the lens and the disk position angle 
are constrained, we find no strong constraints on the disk 
orientation relative to the direction of motion (Figure [5]). 
The convolved uncertainties in the two quantities are too 
large. 

Luminosity that originates outside the accretion disk 
proper is a concern for size estimates because contamina- 
tion by emission on large scales dilutes the microlensing 
signal and le ads us to overest i mate the projected ar ea of 
the disk (see IDai et all [2001 IMorgan et ai. 2008bj). To 
examine this effect, we did a set of trials with (M/M©) = 
0.3 where 0%, 20%, and 40% of the sou rce light was not 
microle nsed. The spectral analysis by lEigcnbrod et al.l 
(2008a) suggests that the level of contamination is 20%. 
The 0% and 20% cases fit equally well, while the 40% case 
had a relative probability 38% lower. As the dilution is 
increased from 0%, more face-on disks are preferred, with 
cosi > 0.63, 0.73, and 0.80 (68% confidence) for 0%, 20%, 
and 40% dilution. Adding unmicrolensed light affects the 
disk size both through the general dilution and the shift 
towards more face-on orientations. The de-projected ra- 
dius is smaller by ~ 26% (55%) if the contamination is 
increased to 20% (40%) from 0%. These corrections are 
not large compared to our statistical uncertainties of or- 
der 70%, but they are an important physical considera- 



tion. 



4. DISCUSSION 

By including random stellar motions in our microlens- 
ing analysis of Q2237 we find evidence that the accre- 
tion disk of Q2237 is viewed face-on with cosi > 0.63. 
This lends support to the popular AGN unification model 
where we expect Type 1 quasars like Q2237 to be viewed 
nearly face-on. Including 20% contamination from broad 
line e mission on larger spatia l scales (estimated from 
spectra lEigenbrod et al.ll2008a[ ) results in a stronger case 
for face-on solutions. Modeling stellar motions and the 
inclination and position angle parameters further reduces 
the systematic uncertainties of our measurements of the 
disk radius compared to earlier studies by including a 
broader range of physical uncertainties. As we found 
in Paper I for the lens velocity and mean stellar mass, 
the results of the separate analyses of LCI and LC2, 
the first and second temporal halves of the OGLE light 
curves, produce consistent results for every parameter 
we considered. While we have used a relatively simple 
model for the accretion disk, these results demonstrate 
that disk shapes can be measured w i th qua sar microlens- 
ing, as suggested by lCongdon et all (|2007| ). As data sets 
and computing power improve, it will be natural to try 
fitting more subtle disk features such as the asymmetries 
from relativist i c effe cts using relativistic models such as 
IHubenv et all (l20?)l . 

The exceptionally long OGLE light curve and the 
fast microlensing timescales of Q2337 made it a natural 
first choice for including the random stellar motions and 
studying the shape of the disk. Our expanded method 
is similar in computational cost to our previous efforts 
with static patterns, so there is no reason not to use it 
generally. The stellar motions clearly aid in reducing un- 
certainties in the mean mass (see Paper I) and allow us to 
correctly make inclination corrections. It can also easily 
be extended to take advantage of multi-wavelength data 
to try to constrain the temperatur e profile. However , 
as in earlier micro l ensing studies (jPooley et all 120071 : 
IMorgan et all [20101 IDai et al.l [2009) , we cannot recon- 
cile the basic temperature profile of a thin disk, the mi- 
crolensing size estimate, and the observed optical flux. 
A disk with a T cx i?~ 3 / 4 temperature profile normalized 
by the microlensing size estimate should be brighter than 
observed. This c an be solved by altering the temper- 
ature profile (see iPoindexter et al.l 120081 : IMorgan et al.1 
2010). For example, reducing the slope from T cx i?r 3 / 4 
to T cx R- 1 / 2 would increase the flux size by a fac- 
tor of 2.2 relative to the half light radius. This would 
be mildly incon sistent (1.5cr) with the sl ope estimate of 
-0.83 ± 0.21 bv lEigenbrod et all (|2008b[ ). However such 
changes also call into question the basic structure of the 
thin disk model. The other simple possibility is to re- 
duce the emissivity of the disk to be well below that of 
a blackbody (by the ratio of the flux/microlensing sizes 
squared), but this seems less physically plausible than 
change in the temperature structure. 

This work was supported in part by an allocation of 
computing time from the Ohio Supercomputer Center. 



7 



REFERENCES 



Agol, E. 1997, Ph.D. Thesis, 

Agol, E., Gogarten, S. M., Gorjian, V., & Kimball, A. 2009, ApJ, 
697, 1010 

Anguita, T., Schmidt, R. W., Turner, E. L., Wambsganss, J., 
Webster, R. L., Loomis, K. A., Long, D., & McMillan, R. 2008, 
A&A, 480, 327 

Antonucci, R. 1993, ARA&A, 31, 473 

Bate, N. F., Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 

2008, MNRAS, 391, 1955 
Blaes, O.M., 2004, in Les Houches Summer School LXXVIII 

(Springer: Berlin) 137 
Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 
Chartas, G., Kochanek, C. S., Dai, X., Poindexter, S., & Garmire, 

G. 2009, ApJ, 693, 174 
Congdon, A. B., Keeton, C. R., & Osmer, S. J. 2007, MNRAS, 

376, 263 

Dai, X., Kochanek, C. S., Chartas, G., Kozlowski, S., Morgan, 
C. W., Garmire, G., & Agol, E. 2009, larXiv:0906.4342l 

Eigenbrod, A., Courbin, F., Sluse, D., Meylan, G., & Agol, E. 
2008a, A&A, 480, 647 

Eigenbrod, A., Courbin, F., Meylan, G., Agol, E., Anguita, T., 
Schmidt, R. W., & Wambsganss, J. 2008, A&A, 490, 933 

Elvis, M. 2000, ApJ, 545, 63 

F loyd, D. J. E., Bate, N. F., & Webster, R. L. 2009, 

[Sxiv:0905.265l] 
Fritz, J., Franceschini, A., & Hatziminaoglou, E. 2006, MNRAS, 

366, 767 

Gil-Merino, R., Wambsganss, J., Goicoechea, L. J., & Lewis, 

G. F. 2005, A&A, 432, 83 
Gould, A. 2000, ApJ, 535, 928 
Hinshaw, G., et al. 2009, ApJS, 180, 225 

Hubeny, I., Blaes, O., Krolik, J. H., & Agol, E. 2001, ApJ, 559, 
680 

Huchra, J., Gorenstein, M., Kent, S., Shapiro, I., Smith, G., 

Horine, E., & Perley, R. 1985, AJ, 90, 691 
Kochanek, C. S., 2004, ApJ, 605, 58 

Kochanek, C. S., Dai, X., Morgan, C, Morgan, N., Poindexter, S. 
& Chartas, G., 2007, in Statistical Challenges in Modern 
Astronomy IV in Statistical Challenges in Modern Astronomy 
IV G. J. Babu and E. D. Feigelson, eds., (Astron. Soc. Pacific: 
San Francisco), astro-ph/0609112 

Kollmeier, J. A., et al. 2006, ApJ, 648, 128 

Koratkar, A & Blaes, O. 1999, PASP, 111, 1 

Morgan, C, Kochanek, C.S., Morgan, N.D., Falco, E.E., 2010, 
ArXiv e-prints, 707, arXiv:0707.0305 in prep. 



Morgan, C. W., Eyler, M. E., Kochanek, C. S., Morgan, N. D., 
Falco, E. E., Vuissoz, C, Courbin, F., & Meylan, G. 2008, ApJ, 
676, 80 

Morgan, C. W., Kochanek, C. S., Dai, X., Morgan, N. D., & 

Falco, E. E. 2008, ApJ, 689, 755 
Mortonson, M. J., Schechter, P. L., & Wambsganss, J. 2005, ApJ, 

628, 594 

Mosquera, A. M., Munoz, J. A., & Mediavilla, E. 2009, ApJ, 691, 
1292 

Peng, C. Y., Impey, C. D., Rix, H.-W., Kochanek, C. S., Keeton, 
C. R., Falco, E. E., Lehar, J., & McLeod, B. A. 2006, ApJ, 649, 
616 

Poindexter, S., Morgan, N., & Kochane k, C. S. 2008, ApJ , 673, 34 
Poindexter, S. & Kochanek, C. S. 2009, arXiv:0910.3213 
Pooley, D., Blackburne, J. A., Rappaport, S., & Schechter, P. L. 

2007, ApJ, 661, 19 

Shakura, N.I. & Sunyaev, R.A., 1973, A&A, 24, 337 
Tinker, J. L. & Wetzel, A. R., Zehavi, I., 2009, in prep. 
Trott, C. M., Treu, T., K oopmans, L. V. E., & Webster, R. L. 

2008, larXiv:0812JJ748l 

Udalski, A., et al. 2006, Acta Astronomica, 56, 293 
Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 
Vakulik, V. G., Schild, R. E., Smirnov, G. V., Dudinov, V. N., & 

Tsvetkova, V. S. 2007, MNRAS, 382, 819 
Vestergaard, M. & Peterson, B.M. 2006, ApJ, 641, 689 
Wambsganss, J, 2006, in Gravitational Lensing: Strong Weak and 

Micro, Saas-Fee Advanced Course 33, G. Meyl an, P. N orth, P. 

Jetzer, eds., (Springer: Berlin) 453, astro-ph/0604278 
Wyithe, J. S. B., Webster, R. L., & Turner, E. L. 1999, MNRAS, 

309, 261 

Wyithe, J. S. B., Webster, R. L., & Turner, E. L. 2000, MNRAS, 
315, 51 

Wyithe, J. S. B., Webster, R. L., & Turner, E. L. 2000, MNRAS, 
318, 762 

Yee, H.K.C. & De Robertis, M.M. 1991, ApJ, 381, 386 
Yonehara, A. 2001, ApJ, 548, L127 

Yoo, J., Kochanek, C. S., Falco, E. E., & McLeod, B. A. 2005, 
ApJ, 626, 51 

Yoo, J., Kochanek, C. S., Falco, E. E., & McLeod, B. A. 2006, 
ApJ, 642, 22 



